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Using the anharmonic oscillator and the He atom as examples, we show that 

the present method provides an efficient scheme in estimating both the ground 

and the excited states. We also discuss the limitations of the present method. 
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I. INTRODUCTION 

In undergraduate quantum mechanics, students are taught generally two types of ap- 
proximation schemes, the perturbation theory and the variational calculation. Each method 
is known to have its own advantages and limitations. The perturbation theory is applicable 
only when the perturbing potential is small, whereas the variational method lacks systematic 
means for improvement after the first approximation. Also, it is rather cumbersome and 
difficult to estimate the excited states in the variational method |T|-^. 

In this paper, we present an approximate scheme which combines the two approxima- 
tion methods. It will be shown that the combined method carries the advantages of both 
methods and, at the same time, overcomes the deficiencies mentioned above. Actually, the 
concept of combining the two approximation methods is not entirely new. For instances, 
this concept has been used in lattice dynamics to obtain a better description of phonons 
in crystals 0. In this approach, perturbation calculations were made on the basis of the 
variational self-consistent harmonic approximations. Also, in relativistic field theories, this 
approach was called the perturbation with variational basis and used to calculate effective 
Gaussian potentials and high order terms P,0|. However, to the authors' knowledge, this 
concept has never been discussed in a plain quantum mechanical language understandable 
to undergraduate or graduate students who do not have knowledge of either lattice dynamics 
or relativistic quantum field theory. Therefore, we believe that developing this combined ap- 
proximation scheme for quantum mechanical problems is not only pedagogically important, 
but also useful for real applications in quantum mechanical problems. 

II. PERTURBATION WITH A VARIATIONAL BASIS 

In a variational calculation, we first choose a trial function \1/(A) as a function of a vari- 
ational parameter A for a given Hamiltonian H. Then, the ground state energy is estimated 
by minimizing the expection value, (\E'a|-^|^a) against A. In the present formalism, we also 



first choose a trial function \1/„(A). However, here, the choice of \l/n(A) is hmited to the cases 
where a parent Hamiltonian Hx^ which carries \l/n(A) as exact eigenfunctions, can be found. 
We will discuss implications of this limitation of the theory below. We just assume here that 
such a parent Hamiltonian can be found for a given problem. Then we rewrite the original 
Hamiltonian, H as follows, 

H = H\ + H — Hx 

= Hx + H' , (1) 

where H' = H — Hx is the perturbing Hamiltonian. Clearly, success of a perturbative 
calculation depends on how small H' can be made. In order to obtain an optimum Hx, we 
first determine A through the condition 

S{^n{X)\H\^n{X))=0. (2) 

Here, we note that this condition is not limited to the ground state, but is valid to all 
excited states also. Therefore, generally A will be dependent on the state number n. Thus 
we express A and Hx as A„ and Hx„ respectively. Using the above notations, we carry out 
the standard perturbation calculation of En- 

En= (^„(A„)|//Aj^n(A„)) 

+ {^n{Xn)\H'\^n{Xn)) 



\{^k{Xn)\H'\^n{Xn))\ 



2 



ky^n En (An) — -C/fc (A„) 

+ ■■■ • (3) 

Here E^\Xn) is the nth energy eigenvalue of Hx„ ; 

HxJ^n{Xn))=El^{Xn)\^n{Xn)) • (4) 

Now we compare the above result with the conventional approximation schemes. First of 
all, for n = 0, the first two terms become 

{MXo)\Hxo + i^'l^o(Ao)) = {MXo)\H\^o{Xo)) , (5) 



which is just the conventional variational result of the ground state energy. The third term 
is the second order perturbation term which provides a systematic improvement over the 
simple variational ground state energy. We note that the perturbative correction term is 
negative for the ground state. Thus, it pushes down the variational values of the ground 
state energy for a better agreement with the true ground state energy. However, it should be 
noted that this correction might overcompensate the true ground state energy as observed 
in the anharmonic oscillator example in the next chapter. The approximate energy which 
combines the variational calculation and the second order perturbation correction no longer 
provides an upper bound to the exact ground state energy. The perturbation expansion of 
Eq.(3) is different from the conventional perturbation theory. The basis functions used in 
the perturbation are those obtained through a variational process and the eigenvalues are 
also optimized values. 

Another important aspect of the present theory is that it can be readily applied in eval- 
uating excited states unlike the conventional variational methods. In the usual variational 
calculation, for the calculation of an excited state energy, it is necessary to construct a trial 
wavefunction which is orthogonal to the ground state wavefunction, thus making the proce- 
dure rather cumbersome and difficult [Q. In the present formalism, it may appear that this 
requirement is automatically satisfied because \I'„(A„) is the nth eigenfunction of H\, which 
is orthogonal to '^^{Xn)- However, this point requires a closer examination. Since A„ is 
determined through minimization as given by Eq.(2), it is different for each n. For example, 
Ao should be used throughout in Eq.(3) for evaluation of the ground state energy, whereas 
Ai obtained from Eq.(2) should be used for evaluation of the first excited state. Therefore, 
the requirement of the orthogonality can be satisfied, but, in an approximate sense, because 
\I^i(Ai) is neither exactly orthogonal to the true ground-state wavefunction nor to \I'o(Ao). 
This point will be made more clear in the next chapters through examples. 



III. ANHARMONIC OSCILLATOR 

In this chapter, we apply the above formahsm to the anharmonic oscillator problem 
which has a x"^ term. The Hamiltonian is given by 

H = -— + ---X + bx , (6) 

where b is positive. The conventional perturbative calculation is possible only when b is 
small [Q. In order to apply the present formalism to this problem, we should first find a 
suitable trial function and the corresponding parent Hamiltonian. A convenient choice is 
given by 

mQ 1 _rnn 2 , , 

^n(A) = (— -)4e 2fi- , (7) 

nil 

for the ground state trial wavefunction. Here fi is a variational parameter. The correspond- 
ing parent Hamiltonian, which carries the above trial wavefunction as the exact ground state 
eigenfunction is given by 

H. = ^ + ^x^ . (8) 

2m 2 ^ ^ 

Using this parent Hamiltonian, we rewrite the original Hamiltonian, 

^ p^ mfi2 ^ ^_^ . (9) 

2m 2 2 ^ ^ ^ ^ 

We denote the nth eigenstate of Hn by \nn). Then the expectation value of {nQ\H\nn) is 
given by 

{nn\ H\nn) 

= {nn\Hn\nn) + {n\H'\n) 

= (n^liif^InQ) {n^\x \n^) + b{nn\x \n^) 

= mn +\)- ^^^^^(2- + 1) + ^(2-^ + 2n + 1) . (10) 



Here, we used the standard quantum mechanical results 0; 

{nn\x^\nn) = {^)\Qn^ + Qn + 3) . (11) 

By taking a variation on {nn\H\nQ) , we obtain a relation which determines Qn, 

Vil-u?^^-— =0 . 12 

m^ in + 1 

This equation is valid for any n. Now it is necessary to evaluate the third term of Eq.(3). 
This second order calculation can also be carried out in a straightforward fashion p[ to yield 

\ hU^ 2 64n^l60n^ -336^3- 664^2 -28n- 24 

\Ml}\m^^l^ (2n + l)2 ' ^ ^ 

Here, we used Eq.(12) to express fi^ — {jj\ in terms of fi^ and h. Collecting terms, we obtain 

1 , bh^ ,,64^5 + 160n^- 336^3-664^2-2871 -24 , , 

• (14) 



Ahrtn^Am^ni' (2n + l)2 

Now we can readily evaluate any energy eigenvalues of the anharmonic oscillator given by 
Eq.(6). First, we obtain Qn from Eq.(12) for a given n and substitute this value to Eq.(14). 
In order to compare the present result with the conventional perturbation and the varia- 
tional calculations, we carry out numerical calculations up to the second order using Eq.(14). 
For this purpose, we choose ^2|L = O.beVA'"^ and carry out calculations for various values 
of b. The results for the ground state energy are shown in Table I. The exact ground state 
energy is obtained numerically using the Runge-Kutta-Fehlberg algorithm 0. The result 
shows that the present method is clearly superior to the conventional perturbation theory 
and the variational calculation. First of all, the present method gives highly accurate values 
in the regime where the conventional perturbation theory is not applicable. This is because 
the perturbational basis has been renormalized through the variational process to yield a 



better convergence. This renormalization effect can be also seen from the changing values of 
^rriQ"^. For large b, ^mQ"^ becomes much larger than |?7ia;^, thus making the perturbation 
Hamiltonian, H' = _ ™(. ~^ ) ^2 _|_ ^^4 correspondingly smaller. Now, in order to examine 
this behavior in detail, we show, in Table II, the results of the first and the second order 
calculations of the approximations for b = 0.05eVA~'^. First, we observe that in the ordi- 
nary perturbation theory, convergence is bad already. The first order result of the present 
approximation scheme is same to the variational calculation as shown in Eq.(3) and (5). We 
can see clearly that the renormalized perturbation calculation yields a good convergence. 

We now extend the calculation to the first excited state energy. The results are shown 
in Table III. For b = O.OSeFA"^, we obtain Ei = 5m2412eV which is 100.02% of the 
true value. Here, the variational parameter ^mQf is 0.990354eyA~^ and clearly different 
from ^m^Q. We note that the agreement for the first excited state is again excellent. As 
mentioned earlier, the value of the variational parameter, Q, is different for each level. 

IV. HELIUM ATOM 

We now apply the present method to evaluate the energy eigenvalues of He atom. Eval- 
uation of the ground state energy using both the standard perturbation method and the 
variational calculation are given in any standard quantum mechanical text books |jl|^ • It is 
known that the standard perturbation calculation does not yield satisfactory result because 
of the large coulomb interaction between electrons. The variational approximation gives a 
better result, which will be utilized below. 

The Hamiltonian for a He atom is given by 

2 2 Vp2 y 2 2 

//=ZL + ^_f^_f^ + ^ , (15) 

2m 2m ri r2 ri2 

where Z is 2. In order to evaluate the energy eigenvalues variationally, we choose a parent 
Hamiltonian 

^,. = rf + ii -£!£!-£!£! , (16) 

2m, 2m ri r2 
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where Z* is the variational parameter [Q. With this choice, the perturbation Hamiltonian 
H' is given by 

and H = Hz* + H' . The ground state eigenfunction for this parent Hamiltonian is given by 

$o(^*) = 0ioo(^*,ri)0ioo(^*,r2)x- , (18) 

where x" = ^(| T )l i ) - I i )l T ))• The expectation value of ($o(^*)|^|'^o(^*)) is readily 
expressed as a function of Z* |[l|], 

{<^o{Z*)\H\<^o{Z*)) = --mc^a'^{4:Z*Z -2Z*'^ --Z*) . (19) 

The optimal value of Z* is obtained by minimizing Eq.(19) and given by Z* = 1.6875. 
Substituting the above results to Eq.(3), we obtain 

E„= - -mc'^a\4Z*Z - 2Z*^ - -Z*) 

^ \{^k{Z*)\H'\MZ*))\' 
+ ^ E^^iZ*) - EtiZ*) 

+ ■■■ (20) 

The first term gives the usual variational result, Eg = —5.6953ryd for the ground state energy 
which is 98.077% of the experimental value of -5.8070 [|l|. In calculating the second order 
contribution, we only have to consider excited state wavefunctions which have antisymmetric 
spin parts, since ^o{Z*) is antisymmetric in spins. The general form of the wavefunction 
with antisymmetric spin parts is given by 

^n'u'm'in,r2) = A[(j)nlm{ri)(pn'l'm'ir2) + (I)n'l'm'iri)(l)nlmir2)]x' , (21) 

where A = 4= if nlm ^ n'l'm' and | if nlm = n'l'm' . For calculation of the second order 
term, we divide the perturbation Hamiltonian, H' , into two parts 

{Z - Z*)e^ {Z-Z*)e^ 



H[ 



n r2 
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and 



e^ 



H' = — . (22) 

ri2 

The H[ contribution does not have any angular dependence. Therefore, it can be straight- 
forwardly calculated to yield 

= —2A{Z — Z*)e ^imfiO^iijn' fioiXn^n' ,1 + -^n'5n,l) ! (23) 

where 

X„= f drr^R^o{r)-R,oir) . 
J r 

The H'2 contribution can be calculated by expanding , ^^ , in terms of spherical harmonics 
0. A lengthy but straightforward calculation yields the result 

i^n'Tm'l H'2 I $100 ) = — 21 + 1 — '^'''"^"^-'"'-^"'^'' ' (^'^) 



where 

ri 
Collecting the above results, we obtain an expression for the ground state energy. 



ynn'i = / drdr2rlrlRni{ri)Rn'i{r2)-^Rio{ri)Rio{r2) ■ (25) 



where 



E^= - -mc^a'^UZ*Z - 2Z*^ - -Z*) 

\/(f)nlm I rr/ 1 ^100\ 12 
I V^ W^nn-ml-"- 1^100/ 1 /r,^\ 

+ 4^ _W2(Z*a)2[2-4-^] ' ^ ^ 

n,n,l,m 2 \ 'I- n^ n i 



l^nlm I TTi I ff,100 \|2 

= I — 2A{Z — Z*)e bimfloiXn^ri' ,1 + Xni5n,l) 

I 2Ae^(-l)- 2 .271 

H oTTl 'Jlm,U-mi'^nn'l\ ■ \^() 



Here, we assume n' > n, < / < n — 1, and < m < I. Also n' = n = 1 should not be 
included in the second order correction. The above expression is numerically evaluated upto 
n' = 7. The second order correction gives —0.02A9ryd, and, thus, makes the total ground 
state energy —5.7202ryd, which is 98.505% of the experimental value of -5.8070. Comparing 
to the original variational calculation value of -5.6953, we observe that a systematic improve- 
ment has been achieved, although the agreement is not as good as the case of anharmonic 
oscillator. We believe that this rather slow convergence may originate from the fact that we 
used the same Z* for both electrons. It is known that if one electron draws closer to the 
nucleus at some instant, it tends to push the other farther out, thus making the screened 
charges different for the two electrons fl]]. Also using the same Z* for excited states may 
introduce additional errors. However, either introducing an additional variational parameter 
or calculating Z* for each excited level would make the present scheme too complicated to 
be useful. We just note that for the first excited state of symmetric spin state, Z* is given by 
1.8497 which is quite different from that of the ground state. The first excited state energy 
obtained variationally without the second order correction is given -4.2765 ryd comparing 
to the experimental value of -4.3504 ryd. Because of computational complexity, we have 
not carried out the full second order calculation. However, partial calculations with some 
selected terms also show the same trends as in the case of the ground state. 

V. CONCLUSIONS 

We have presented an approximation scheme which combines the perturbation theory 
with the variational calculation in quantum mechanics. It is shown that this scheme provides 
a very good convergence beyond the first variational calculation even when the perturbing 
potential is large. Also it can be readily used to estimate excited states. Therefore, this 
method overcomes the shortcomings of the perturbation theory and the variational calcu- 
lation and combines the advantages, when applicable. Here, it should be noted that the 
present method is applicable only when a parent Hamiltonian can be found, thus limiting 
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its usefulness. However, it is known that for a quite large number of problems, potentials 
can be expanded into simple harmonic potential plus higher order terms. In such category 
of problems, the present method is expected to be efficient. 

Another merit of the present method which, we believe, may be more important is the 
pedagogical value of the combined approximation scheme. The theory shows to the students 
the limits of both approximations and demonstrates a way to overcome those, although in 
a limited class of problems. Also, this method provides an easy example of renormalized 
perturbation theory which are often used in many body and field theories JlO . 



ACKNOWLEDGMENTS 

This work has been partially supported by the Korean Ministry of Education (BSRI-97- 
2425), the Korea Science and Engineering Foundation through Project No. 95-0701-04-01-3 
and also through the SRC program of SNU-CTP, and Yonsei University Research Fund. 



11 



REFERENCES 

[1] A. Goswami, Quantum Mechanics (Dubuque: Win. C Brown, 1992) pp 388-395. 

[2] R. L. Liboff, Introductory Quantum Mechanics (New York: Addison- Wesley, 1997) 3rd 
ed. pp 640. 

[3] S. Gasiorowicz, Quantum Physics (New York: Wiley, 1995) 2nd ed. pp 296-304. 

[4] R. Shankar, Principles of Quantum Mechanics (New York: Plenum Press, 1994) 2nd 
ed. pp 429-435. 

[5] T. R. Koehler, Phys. Rev. 165, 942 (1968). 

[6] P. Cea and L. Tedesco, Phys. Rev. D 55, 4967 (1997). 

[7] G. H. Lee and J. H. Yee, Phys. Rev. D 56, 1 (1997). 

[8] S. Fliigge, Practical Quantum Mechanics (Berlin: Springer- Verlag, 1971) pp 80-84. 

[9] P. L. DeVries, A First Course in Computational Physics (New York: Wiley, 1994) pp 
219. 

[10] J. W. Negele and H. Orland, Quantum Many-Particle Systems (New York: Addison- 
Wesley, 1988) pp 373. 



12 



TABLES 
TABLE L Comparison of the ground state energies (in eV) obtained using different approxima- 
tion schemes. The values in the parentheses are the ratios to the exact values. The values of ^^^-^o 
are also shown as references. 



b(eVl-4) 


0.01 


0.05 


0.25 


Perturbation 
Theory 


1.4318427 
(99.935%) 


1.5279252 
(95.962%) 


does not converge. 


Variational 
Calculation 


1.4333279 

(100.038%) 


1.5968858 
(100.293%) 


2.0664772 

(100.929%) 


Present 
Method 


1.4327276 

(99.997%) 


1.5912088 
(99.937%) 


2.0412648 
(99.679%) 


Exact 
Energy Value 


1.4327725 


1.5922195 


2.0474629 


^mnlieVA-^) 


0.5770839 


0.8227827 


1.6423320 



TABLE IL Detailed comparison of the ground state energy values (in eV) obtained by different 
methods for b=0.05 eVA~^. 



Method of 
Approximation 


First order 
calculation 


Second order 
calculation 


Perturbation 


1.6659633 
(104.632%) 


1.5279252 

(95.962%) 


Variational 
Calculation 


1.5968858 
(100.293%) 




Present 

Method 


1.5968858 
(100.293%) 


1.5912088 
(99.937%) 
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TABLE III. Comparison of the first excited state energies (in eV) for b=0.05eVA . The values 
in the parentheses are ratios to the exact value. The value of 2^1^! is shown together as a reference. 



Perturbation 


4.484801(88.09%) 


Variational Calculation 


5.106102(100.29%) 


Present Method 


5.092412(100.02%) 


Exact Value 


5.091282 


iml7f(eFl-2) 


0.990354 
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